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Abstract 

A method is proposed to construct spiral curves by inversion of a spiral arc of parabola. The 
resulting curve is rational of 4-th order. Proper selection of the parabolic arc and parame- 
ters of inversion allows to match a wide range of boundary conditions, namely, tangents and 
curvatures at the endpoints, including those, assuming inflection. 
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1. Introduction 

Constructing curves with monotone curvature function (spirals) is a well-known problem 
in CAD applications. One application of finding spiral arcs with predefined boundary condi- 
tions is the design of transitions curves, joining two given curves in G^-continuous manner. 

Another application was proposed in [11]. For a spiral, represented by a set of interpolation 
nodes with boundary tangents, a region, covering all possible instances of spiral splines, was 
constructed. The width of the region allows to estimate the determiness of a curve by 
the given point set, and allows to accord curve design with the requirements, imposed by 
manufacturing tolerances. 

The recent approaches to the problem include approximate solutions for small arcs, solu- 
tions for specific boundary conditions, attempts to select spiral segments within traditional 
polynomial or rational curves. An exact solution was given in ^ in terms of Cornu spiral, 
extended by two circular arcs. 

The author's contribution to the subject is the study of spirals as a special class of planar 
curves jsl, Q, [sj . The main results of this study are the necessary and sufficient conditions 
of existence of a spiral with given curvature elements at the endpoints. They expand some 
statemets of W. Vogt and A. Ostrowski, well known from |6j], onto non-conex spiral arcs, i.e. 
those with inflection or multiple winding. The requirement of curvature continuity was also 
removed, which is the case of circular splines with piecewise constant curvature. 

Based on this study, a construction in terms of rational curves was developed and is being 
presented in this note. The construction is straightforward and does not require any heuristic 
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optimizational or fitting procedures. A wide range of boundary conditions is covered. A 
possible development of tfie metliod to satisfy any boundary conditions, compatible with 
spirality, is discussed in the last section. 

1.1. Overview of the method 

The proposed method is illustrated by Fig. [H Suppose, we have to construct a spiral 
arc from the point Ai = {—1,0) to i?i = (1, 0). Required parameters are shown as tangent 
vectors and the circles of curvature (dashed) at the startpoint Ai, marked by a and k^, and, 
at the endpoint Bi, hj j3 and k^. 

First, we construct a specially calculated parabolic arc, shown by dotted line, inscribed 
into its control polygon AiPiBi. This is a second order Bezier curve 

{x{t),y{t))=Ai{l-tf + 2Pi{l-t)t + Bit\ O^t^l, (1) 

with the control point Pi = {p, q). At first glance, this arc has nothing in common with the 
required curve. The only visible feature is that it is a spiral: the vertex of parabola, although 
close to the point Ai, is outside the arc AiBi, t^ ^ (0, 1). 

Second, we apply linear fractional transformation (also known as Mobius transformation) 

Zo) = T~r~^^ ^0 = XQ+iy^ 7^ ±1, (2) 

i + ZqZ 

where z = x{t)+iy{t) is a point on the parabola, Zq is a complex constant. The image of 
the parabolic arc is the sought for curve. In this example we have calculated p ~ —0.8845, 
q ~ -0.3033, xq ~ 1.0296, y^ ~ -0.6727 for prescribed values a = -180°, k^ = 2.5, (3 = 120° 
and kg = 0.5. 

The solution is based on the following features: 

• we can vary 4 free parameters, {p, q, xq, yo}, to satisfy 4 required values {a, (3, k^, k^}; 

• Mobius transformation preserves monotonicity of curvature; if {p, q) defines a spiral 
parabolic arc (without the vertex inside), the transformed curve is a spiral; 
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• original curve ([T]) being 2nd order polynomial, the image is 4th order rational. 

It is well known [3] that transformation ([2]) includes movements, homothety, symmetry 
and inversion. The latter provides the necessary flexibility in modifying the form of a curve. 
In Fig. [T]& a curve with inflection was requested {a = (3 = —40°, 3 = k^>0> kg = —2) and 
constructed by similar deformation of the parabola with the control point P2 . 

2. Definitions and notation 

An arc of a curve is described by the functions of the parameter t: 

x{t), y{t), z{t)=x{t)+iy{t), T{t), kit), O^t^l, 

r and k being the angle of the tangent vector and the curvature at the point (x, y) . In this 
article we deal with curves whose function k{t) is monotone on the segment t G [0, 1], i.e. 
with spiral arcs. 

Definition 1. A spiral arc AB is short, if its tangent vector never achieves the direction 
BA, opposite to the direction of its chord, except, possibly, at the endpoints. 

Another definition, namely, that "a spiral arc is short, if it does not intersect the com- 
plement of the chord to the infinite straight line (possibly intersecting the chord itself)", is 
equivalent to Def.[T] 0, lemma 4]. Four spiral arcs AiMiBi, z = 1, 2, 3, 4, in Fig. [2]are short. In 
this article we consider short spirals only, which is sufficient for most of CAD applications. 
Examples of long ones are: 
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We denote below a circle of curvature as a quadruple of coordinates of a fixed point, 
tangent and signed curvature at this point: 

It may be a straight line, if ki = 0. 

To consider properties of the arc AB with respect to the chord AB of the length \AB\ = 2c, 
we choose the coordinate system such that the chord becomes the segment [— c, c] of the x- 
axis. With a = r(0), k^ = k{0), and /5 = r(l), /c^ = A;(l), the boundary circles of curvature 
take form 

/Ci = /r(-c,0,a,fcJ, }C2 = K{c,0,p,k^). (3) 

It is convenient to assume homothety with the coefficient c~^, and to operate on the segment 
[—1,1] of the X-axis. The coordinates x, y, arc length s, and curvature k become normalized 
dimensionless quantities, corresponding to x/c, y/c, s/c and kc. With such homothety 
applied, boundary circles ([3]) appear as (jl]). 

Definition 2. A spiral arc with boundary curvature elements 

/Ci = /^(— 1, 0, a, a), /C2 = -ft'(l, 0, /3, 6), a = ck^, b = ck^. (4) 

is said to be in normalized position. The product kc, invariant under homotheties, will be 
referred to as normalized curvature. 

For each curve AiMiBi in Fig. [2] two circular arcs are traced from Ai to Bi. They share 
tangent with the spiral at one of the endpoints. These two arcs form a lense, enclosing the 
spiral. The angular width of the lense (signed) is 

a = a + j3. 

Four points Mj in Fig. [21 four arcs A^MiBi (the first one is parabola) and four lenses are 
images of each other under transformation ([2]). The angle a remains constant for the four 
lenses. 

3. Theoretical background 

The following is valid for spirals: 

(i) Monotonicity of curvature is preserved under inversion, increasing curvature being 
transformed into decreasing, and vice versa j^, theorem 1]. 

(ii) The necessary and sufficient condition for the existence of a non-biarc spiral with 
boundary conditions ([3]) is inequality 

Q(/Ci, IC2) = {k^c + sin a){kgC — sin 13) + sin^ — < (5) 
, th. 2]. If Q = 0, the biarc is the unique spiral, matching these boundary parameters. 
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(iii) For the existence of a short spiral, it is required additionally that 

sign(a + /3) = sign(A;^ - k^) ^ 0, 

0, th. 10, or, in details, 

if < kj^ : a+P > 0, — vr < a, P ^ n; , , 

if k^> kg : a+[3 < 0, -vr < < vr. ^ ' 

(iv) Both Q and a = a+(3 are invariant under Mobius transformations. 

(v) Short spiral is enclosed in its lense [sl, theorem 2]. 

Below we provide some comments to these statements. 

Alternative proof of ^ could be easily derived from the well known facts that (a) the curve 
in the vicinity of some point intersects its circle of curvature if the curvature is monotone, and 
(b) the curve is located from one side of this circle, if the point corresponds to vertex. These 
local properties are evidently invariant under inversions, and no new vertices can appear on 
the transformed curve. 

In (In]) and ^ Q = sin^ ^, where is the intersection angle of the two circles of curvature. 
Q < means that ip is pure imaginary. Q < can be interpreted as follows: if two given circles 
are inverted into a concentric pair, taking into account their orientation, the resulting circles 
will be parallel (and not anizparallel, as it will be the case ii Q>1). Q = Q involves ^ = 0, 
and corresponds to tangency of two circles (biarc curve can be constructed). < Q < 1 means 
their real intersection, which is not compatible with spirality. Neither is Q ^ 1- 

Statement dm]) is illustrated by Fig. [1] for decreasing curvature (a = a+[3 < 0). Note that 
if we bring the function r(t) to the range [— tt, tt], then, to preserve continuity, the value — tt, 
not +7r, should be assigned to the angle a of the spiral AiBi in Fig. [T]a. This accords to (JH]). 
Curvatures k{t) in Fig. [2] are increasing, and a > 0. 

For (jivT) we note that oriented angles if) and a only change sign under inversions, where- 
from the invariance of Q, \a\ follows. Transformation ([2]) is either identity (zq = 0), or includes 
both inversion and symmetry. The type of monotonicity of curvature (increasing/ decreasing) 
is swapped under inversion and restored under symmetry. So does sign a. 

Statement (jvj) was proven for convex spirals in [sl and for any short spiral in j^. Fig. [2] 
demonstrates it, including non-convex case A2M2B2. 

Fig. [3] illustrates inequahty Q < ([5]) , for fixed angles a, (3 taking form 

Q(a, b]a, P) = {a + sin a) (6 — sin /3) + sin^ + P ^ 

Angles a = 45°, /5 = -15° are those of the arc AiMiEi, the point 7^^= (a, 6) ~ (-1.98, -0.10) 
corresponds to normalized boundary curvatures of this arc. The equation Q{a,b) = de- 
scribes a hyperbola in the plane (a, b). Two convex regions Q ^0, bounded by two branches 



-'^ Articles [1, are not yet translated in English. The interested reader may find English version of all 
these facts and proofs in arXiv: math/060144C)| ^2 (A. Kurnosenko. Around Vogt's theorem). 
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of the hyperbola, cover the values of boundary curvatures, allowing spiral arcs to be con- 
structed. The upper left branch, lying in the halfplane a<b, corresponds to increasing 
curvatures. As the angles in the first example are such that a+(3 = 30° > 0, corresponding 
short spirals are of increasing curvature IQ; therefore, possible boundary curvatures (a, 6) 
of short spirals for this case are covered by the upper left region. 

The second example in Fig. [3] is drawn for a = —180°, (3= 120°, the angles of the spiral 
AiBi in Fig. [T]a. As cr = —60° < 0, short curves are of decreasing curvature, a>b. Possible 
values of (a, b) are covered by the lower right hyperbolic region. 

The following inequalities just reflect the position of two regions in question with respect 
to the asimptotae of the hyperbola Q, cor. 2.1]: 

a^b =^ a + sina^O, 6 — sin /3 ^0. (7) 

Curvatures at the center of the hyperbola (ao = — sin a, bo = sin P) are those of two circular 
arcs, bounding the lense. 

Additional hyperbola, drawn dashed, bounds the region of applicability of the proposed 
method. This is discussed later (Prop.H]). 

4. Transformation of a spiral arc 

Now we investigate map ([2]). This is a particular case of general Mobius transformation 
z — >• ^f^, ad— be 7^ 0, keeping points (—1, 0) and (1, 0) intact. 

Proposition 1. Let two eircles of curvature @ of a spiral are, and two eireles 

K.\ = K{-1, 0, a\ a""), K.\ = K{1, 0, (3\ 6*), 
of another spiral are are sueh that 

cr(/Ci,/C2) = (T(/ct,/c^), g(/Ci,/C2) = g(/ct,/c5). (8) 
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Then transformation ([2]) with 

r Qi\o _ 1 f Ao = a*-a = /3-/3* , 

zo = — , where < a + sin a 6*-sin/3* (9) 

a* + sm a o — sm p 

ma;?5 the fist pair of circles to the second one. 
Proof. Conditions (IHD can be written as 



a + p = a* + p*, 

(a + sin a) (6 — sin /3) = (a* + sin a*) {b* — sin /3* 



(10) 



that's why two expressions for Aq and ro in are equivalent. Because sign(a+/3) = 
sign(a*+/?*), two spirals have the same type of monotonicity of curvature. Numerator and 
denominator in ro are non-zero and of equal signs ([7]). Thus ro > in ([9]) is justified. 

Treating z in W{z; zq) as any object, subjected to map ([2]), we denote images of circles 
/Ci,2 as W{]Ci^2', zq)- We have to prove that W{]Ci^2', zq) = A^i 2- 

Let z{s) be arc-length parametrization of the circle /Ci, i.e. 

z(s) = -1 + -e^"(l - e''^^) : z',{s) = e'("+'^^), = iae'("+'^^). 

a 

We need not pay special attention to the straight line case a = 0, because 

lim z{s) = -1 + se'" 

is the straight line in question. Let w{s) = W{z{s); zo) be some parametrization of the image 
of /Ci (in general case the parameter s is not arc length for w{s)). Calculating derivatives at 
z{0)=w{0) = -l yields: 

w'As) = zi{s)- 



[l + zozis)]'' 

<(0) = ^s(0) n i~?n^i2 = e^"^^ = ■ roe'^" = ro&^\ 
[1 + zoz{0)\^ I- Zo 

As ro > 0, we get required tangent a* for w{0). To control curvature, we need the second 
derivative: 



[l + z,z{s)f ' ' [l + Zoz{s)]^' 
w'^^{0) = iae*"^^^ - 2^06^'" ^ \ = roe'°*(ia + e^" - roe'"*). 



We separate real and imaginary parts, w^(0) = ui + ivi, Wgg{0) = U2 + iv2: 

ui = ro cos a*, vi = ro sin a*, 

U2 = ro [cos(a+a*) — a sin a* — ro cos(2a*)] , 

V2 = — ro[sin(Q;+a*) + a cos a* + rosin(2a*)]. 
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Calculating the curvature of the circle w{s) at s = proves that W{]Ci; Zq) = /C*: 

V2U1 — V1U2 a + sin a 



— sm a = a 



Similar calculations for the second pair of circles can be omitted. The required values 
(3* and h* result from ([8]) and from the invariance oi Q,a (jiv|) . Namely, denote the image 
W{1C2]Zq) as if(l,0,/?,6). Then 



Thus obtained equation a (/C^, W{1C2] zq)) = o"(/C^, /C2) is in fact a*+/3 = a*+P*, and yields 
P = P*. The same reasoning for Q yields b = b*, i.e. W{)C2; zq) = IC^- □ 

5. Spiral parabolic arc 

In this section we explore a spiral parabolic arc from the viewpoint of its boundary 
conditions and invariants. The equation ([1]) with the control point (j9, q) can be rewritten as 

x{t) = -{l-tf + 2p{l-t)t + t\ 

y{t) = 2qt{l-t), U^^^^- ^'') 



Calculate 



x'{t) = 2 (l+p)(l-t) + 2{l-p )t, y'{t) = 2q{l-t) - 2qt, 

g{t) = ^x''\t)+y''\t) = 2^{l~t)^hl + {l-t)t{hl+hl-A) + t^hl 



where 



= \AP\ = v/(l+p)2 + g2 and h2 = \PB\ = ^(1-^)2 + ^2. (12) 
Since cosr(t) = x'(t)/g{t), sinr(t) = y'(t)/g(t), the boundary angles are defined by 

l+P . q n '^-P ■ n -Q /I ox 

cosa = — — , sma; = — , cosp = — — , smp = - — , (13) 

hi hi h2 Ai2 



and invariant a = a+f3 by 



l-p'^ + q'^ . -2pq 

cosa = ; — ; , smcr=- — ; — . (14) 

hih2 hih2 



Curvatures are 



I f+\ y'x'~xy -8q -g -g 

Mt) = ^^ = ^, « = M0) = ^, b = kil) = ^. (15) 

Invariant ([5]), expressed as a function of the control point, looks like 

n( \ 1 , p' + q'-l q\2p^ + 2q' + l) 

^^^'^) = 2+^M^ hm — ■ ^ ^ 
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Figure 4. 



Proposition 2. Parabolic arc (ITTl) is spiral if and only if the control point {p, q) satisfies 
inequalities 

F,{p,q)-F2ip,q)^0, q ^ 0, 
where Fi{x,y) = x"^ + x + y'^, F2{x,y) = x'^ — x + y'^. 

Proof. Differentiating of k{t) f|T5|) yields 

I92q\(p'^ + q^)(2t-l) -p] 192g , , 
k'{t) = '^^^ '—^ = ^ [-(1 - t)F,{p, q) + tF^ip, q)] , 

k'{0) = -^F,{p,q), k'{l) = ^F2{p,q). 

To avoid vertex to occur in t G (0, 1), we require tliat tlie derivative did not cliange sign 
witfiin tliis interval, allowing the vertex at t = or t = l: k'{0)k'{l) '^0. This yields the 
first inequality in f|T7j) : the second one prevents parabola from degenerating into a straight 
line. □ 

Fig. Hlfl illustrates (|T7j) : the control point should be taken within or on the boundary of 
any of two circles Fi^2{x,y) = 0, but not on the x-axis. 

Proposition 3. The locus of control points {p, q), yielding o = const, is 

H{x, y) = sin o"(l — + y^) + 2xy cos cr = 0, sign{xy) = — sign a. (18) 

Proof. The proof follows immediately from (I14p . □ 

The locus (fTSi) is an equilateral hyperbola, centered at the origin, passing through points 
A and B (Fig. |1]&). The second equality in (fT8|) keeps only the part of the hyperbola, lying, 
if (7>0, in quadrants II, IV. Finally, subarcs AAi and BiB provide control points {p,q), 
corresponding to spiral parabolic arcs. For cr < we get the picture, symmetric about the 
X-axis, with subarcs AAi and BiB in quadrants I, III (Fig. |l]c). 
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Both conditions (fTSjl in polar coordinates (p, ^ook like the polar equation of the hyper- 
bola and the intervals for ^: 

ysin(a-20 ^<0: o<e<- + -, -7r<e<---. 
Proposition 4. Under conditions 

0<|ao|<|, (20a) 



Qo ^ Qma. = "^'^f^VK^. ^ ^ = f/tan^. (20b) 



w;2)(w;2 + l)3 

t/iere exist two parabolic spiral arcs fllip . smc/i t/iat Q(j9, (?) =Qo c^'^'^ cr(p, Q') =o'o. 

Proof. The direction of tangent to the hyperbola H{x,y)=0 at points A and B is ctq 
(Fig. |1]6). Condition fl20ap assures the existence of the non- vanishing arcs AAi and BBi 
within limiting circles f[T7|) . 

Consider the behavior of the invariant Q{p, q) while the control point moves along the 
path AAi (or BBi). Q{p,q) can be considered as the function Q{^) of the polar angle. To 
get it, we first express the product /11/12 from (fT4|) : 

—2pq — 2p^(^) sin,^ cos,^ sin2,^ 



sin (To sin (Jo sin ( 2^ — ctq ) 

Substituting p^ + q^ = p^^C), q^ = P^(0 sin^ and the above product into f|T6l) . we obtain 

Q(^; o-q) =/2(o^o) - sin^ o-o/i(0> 

3 1 

where /2(a) = cos^ (T cosaH — , /r,-,^ 

2 2 l^-Lj 

, , l + 4cos2^-8cos^{ tan^^ + etan^^-S 

and Ji{q,) = 5 = ■ 

8 cos"* ^ sin ^ 8 tan ^ 

Q(^; (Jq) is a monotone function of 

dQ —3 sin^ (Tq 
d^ ~ 8 cos^ ^ sin^ ^ 

As ^ decreases from A to Ai and farther to the asymptota of the hyperbola, Q increases 
from lim Q(^;(Jo) = — 00 to lim (5(^;cro) = l, 

remaining still negative at Ai, because parabolic arc is still spiral (jSD- We denote this value as 
Qmax- To calculate it, find points Ai, Bi by intersecting the pair of circles Fi{p, q)F2{p, q) = 
with the hyperbola, i.e. from the equations 

{p"^ + q^y - p^ = 0, p = p(Ocos^, g = p(Osin^. 
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They are simplified to p^(0 = cos^^ and, finally, to 

tan f ( tan'^ f tan — + 3 tan £ tan — 
^ V ^ 2 2 



1— tan 



2 c^o 



We ignore the root tan^ = (intersection points A,B), and substitute tan^ = w^: 



tan^ ^ + 3w^ tan ^ + 1 — = 
(lutan,^ + l—w"^) [w^tan^,^ + tu(w^ — 1) tan,^ + 



0. 



The unique real root tan^i 
the expression ( I20b[) for Q„ 



gives both points Ai and Bi. Calculating Q{^i; o"o) yields 



□ 



The plot Qmax{<^) is shown in Fig. [5l The region of curvatures, 
rejected by condition (120bp . was shown in Fig. [3] as the band-like 
zone between the line Q{a, b) = and the dashed line Q{a, h) = Qmax- 
The narrower is the lense, the narrower is the rejection region 
Qmax <Q ^0. In other words, our method is not applicable, when 
boundary circles of curvature are close to tangency. 

Now we find the explicit solution ^ = .^o for the equation 
QiCl'^o) = Qo, i-e. find the control points on hyperbola (fT9|) with 
(T = (Jo, yielding the parabolic arc with predefined value of invariants 

c"o, Qo- 

We denote for brevity 

Qo - /2(o^o) 




Figure 5. 



Qi 



sm (To 



and 



ri = 
r2 = 
ri2 



\l'm—\ 



m 



\Qi 



i + Qi 



n 



\/m^+m+l , 



(22) 



n 



^/2n - (m+2) 



V2n + m+2 



2n - (2m+l) 



2n + 2m + 1 

Each of the above definitions for r's is represented in the alternative form, intended to avoid 
loss of precision while subtracting close positive numbers in calculations. For the same 
reason, the expression for 6^0 is splitted in ( 123!) . and r2—r\ is replaced by ^ 



ri2 



Vr2 



Proposition 5. Solutions ^o,^o of the equation Q{^]<yo) =Qo, o,^d corresponding control 
points (pi,2, 91,2) cire given by 

ri2 

sign (To, if aoQi < 0; 



6*0 = -ri sign Qi - r2 sign (Tq 
^0 = -^o+TT, (P2, ^2) = {-Pi, -gi 



ri + r2 

-(ri + ra) signao, if a^Qi ^ 0; 
^0 = arctan6'o, (pi, qi) = (p(^o) cos^o, p(6) sin^o) , 



(23) 
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Proof. Equation (121]) . namely Qq = /2(o"o) — sin^cro/i(05 can be transformed to 

+ 69^ + 8Qi9 -3 = 0, with 9 = tan ^, 



(24) 



and Qi, defined by (l22l) . Following Descartes-Euler's method, find the cubic resolvent of 
( I24|) and its roots: 



(z + - 64(1 + Ql] 



0, 



^1 

2^2,3 



-4 + 4m, 
-4 + 4me=^i' 



Now one should choose signs for Q = -^^JYi to satisfy C1C2C3 = — SQi. We get it, assuming 

Ci = -signQiyir= -2risignQ i, 
(2,3 = +v/^ = r2 ± iV2n + m+2. 



thus obtaining C1C2C3 = — SsignQiV^'^ — 1 = - 
of (IMD are equal to 

^l,2 = ^Cl±^(C2+C3) 

and are of opposite signs: 

^1^2 = r\- rl 

From (IT^ we note that 

2p^(0sinCcos^ 



sm (Jo 



IsignQi ■ IQil = — 8Q1. So, both real roots 
-ri signQi ± r2. 

-ri2 < 0. 

sign (Jo = — sign tan ,^ = — sign 9. 



So, the only admissible root 9{^ from 6'i^2 is of the sign, opposite to that of (Jq: 



if (To > 0, Qi < : 00 



ri — r2 < 



if (To > 0, Qi ^ 
if (Jo < 0, Qi ^ 
if (To < 0, Qi > 



6'o = -Ti - r2 < 
6)0 = n + r2 > 
6*0 = + r2 > 0. 



This selection is unified in fl23l) . 

Thus found ^ = arctan 6*0 provides two symmetric parabolas, shown in Fig. Hie. 



□ 



6. Algorithm 

Now we summarize step by step the construction of a short spiral arc with predefined 
boundary conditions. 

Step 1. Transformation to the chord's coordinate system. 

Given boundary conditions Km = K{xi,yi,Ti, ki) at the startpoint M (Fig. [6]a), and 
)Cn = K{x2, 1/2, T2, ^2) at the endpoint A^, transform them by orthogonal transformation to 

/Ci = K{-1, 0, a\ a'^), /C2 = K{1, 0, (3\ b*), 
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Figure 6. 
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where 



a* = Ti — jj,, a* = kic 

(3* = T2- fi, h* = k2C. 



/i = arg [(X2-X1) + i(i/2-l/i)] , 



The transformed configuration is shown in Fig. [6f>. 
Step 2. Check solvabihty of the problem. 

• Calculate Qq = QiJCl.lC^) and check inequality ([5]). If (5o>0, boundary conditions 
are invalid, irrespective of the proposed method: no such spiral exists. If Qq = 0, biarc 
is unique spiral. Continue, if Qo < 0. 

• As specified in ([6]), bring the boundary angles to the halfintervals {—tt,tt] or [—71,11). 
Namely, if |a*| = vr or = vr, replace +11 by — vr for the case of decreasing curvature. 
Calculate ctq = check condition ([6]). If it fails, no short spiral exists. 

• Verify conditions fl20|) . If Qq > Qmax or |cro| ^ 7r/2, the proposed algorithm is not 
applicable. 

Step 3. Find spiral parabolic arc ffTTj) such that its boundary tangents a, (3 and curvatures 
a,b satisfy equations ffTU]) . Two solutions, (^1,2, ^'1,2)5 are explicitly given above fl2^ . The 
explicit expression for Qi in ( l22i) is 



Two parabolic arcs are shown in Fig. [6]c,d by dotted lines. 

Step 4. Apply this step to each of two just constructed parabolic arcs. 

Calculate boundary tangents a, (3 from (fT3!) and curvatures a, h from (fT5|) . Define ro, Aq 
from ([H]). Transforming the parabolic arcs by ([2]), get the sought for rational curve: 



where x(t),y(t) are abbreviated to x,y, and If = + ^2 = (x— 1)^ + The second 

expression is tolerant to the case Zo = oo, which is described by finite values tq, Aq. 
Two solutions are shown combined in Fig. [6]e. 

Step 5, final. Return to the original coordinate system (Fig. [6]f). 



Ql = cot (To + 



{a* + sin a*) (6* — sin 
sin^ (To 



X{t) + iY{t) 



1 + 2xxo - 2yyo + {x^ + y^){xl + y^) 



TqII — + 2iro[2?/cos Aq — (x^ + y"^ — 1) sin Aq] 
r2/2 _ 2ro[2y sin Aq + (a;^ + y'^-2) cos Aq] + ll ' 
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Figure 7. 




7. Illustrations 

In Fig. [6] transition curves, joining two concentric circles, have been constructed. Fig. [7] 
adds the curvature plots with respect to arc length for both solutions. 

In Fig. [8] the boundary circles are also concentric, but the tangents at points A and B 
are such that no short joining spiral exist: condition ([6]) is violated. To construct a long one, 
we introduce the point M at the polar halfway from A to B, and new intermediate circle of 
curvature, concentric to both given ones. Two shorts spirals are constructed on two chords, 
AM and MB. Combining them together, we obtain four variants of the long spiral AMB, 



15 





k 

-1 ^-^"^^ 


















1 C, Ci Cs C4 C5 Cfi C7 ^ 

















[q- 
















c^ 










/•/ / 










Cm 






1 • ' — * 

Co 


1 • * — — • 

Co 



Figure 10. 



and the curvature remains continuous at M. In particular, taking radius R^^ = ^JR^R^, 
we get similarity of boundary conditions. Arcs MB can be obtained just by rotation and 
homothety, applied to both arcs AM. 

In Fig. [9] symmetric boundary conditions are given: tangents are parallel, and A;^ < < 
~ I^aI- Curvature changes sign, inflection point on the transition curve is therefore 
required. It does appear in both solutions. 
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Figure 11. 



One could expect a symmetric solution, but, instead, we get the symmetry among two 
solutions. The problem could be resolved, as in the previous example, by introducing the 
intermediate inflection point at the coordinate origin. 

In Fig. [in] the long arc CqCj of Cornu spiral is drawn by dots and subdivided into 7 
short subarcs, CjCj+i. The boundary conditions {xi,yi,Ti, ki) are borrowed from those of 
Cornu spiral at points Cj. In this example every point Ci is taken almost as far from Cj_i as 
limitation ( I20b[) allows. As before, we get two solutions on every segment CjCi+i. Differences 
between two solutions and Cornu spiral itself are well visible on the segment CqCi. This can 
be avoided, as shown on the right side, by inserting an intermediate point C1/2 between Co 
and Ci. 

The curvature plots with respect to arc length for approximations of Cornu spiral are also 
shown in Fig. [10] for both solutions on every segment of the curve CqC-j. Abscissas of vertical 
lines correspond to arc length of Cornu spiral in points Cj. Small vertical marks separate 
arcs of approximating curves. Proportionality k{s) ~ s of Cornu spiral is well reproduced in 
approximat ions . 

8. Conclusions 

Note that we get much more flexibilty in form control by involving spiral arcs of other 
conies or other spirals: basic Proposition [T] is independent of the kind of spiral involved. 
Arcs within a quarter of an ellipse are spirals, and, subjected to transformation ([2]), provide 
another variety of spiral arcs. Being 2-nd order rationals, they also yield 4-th order rational 
spirals. So do hyperbolic arcs. 

Fig. [m shows inversions of a half-hyperbola, traced from the vertex at point A through 
infinity to the second vertex at B. Several inversions of this branch of hyperbola are shown, 
every example with the circle of inversion. Inversion with respect to the unit circle produces 
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a curve AOB, known as hyperbolic lemniscate. Although the original curve has discontinuity 
at infinity, this feature disappears under inversion. The infinite point goes to the center of 
inversion and becomes an ordinary, infinitely differentiable, point of the curve-image. 

Varying eccentricity of conic gives additional possibilities in constructing rational spirals. 
A family of curves is expected as a solution. Starting with hyperbola, one could get a desired 
symmetric solution for the symmetric problem in Fig. [9l 

Preliminary investigation shows that the region of boundary conditions fl2UI) can be ex- 
tended to I (To I ^ vr and Qo<0. Boundary angles of the lemniscale AOB are a = P = n/2, its 
lense is the unit circle, and the width of the lense is a = 11. Taking eccentricity e ^ 1, two 
branches can be made anyhow close to parallel lines (Fig. UTh). The image can be anyhow 
close to a biarc curve, and Q close to 0. 

In this article we have considered parabolic arcs only, thus illustrating the simplest, but 
sufficiently powerful, version of the proposed method. 
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